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Abstract: 

The Square Kilometre Array (SKA) provides an excellent opportunity for low cost searches for fast 
radio transients. The increased sensitivity and field of view of the SKA compared with other radio 
telescopes will make it an ideal instrument to search for impulsive emission from high energy density 
events. We present a high-level search 'use case' and propose event rate per unit cost as a figure of 
merit to compare transient survey strategies for radio telescope arrays; we use event rate per beam 
formed and searched as a first-order approximation of this measure. Key results are that incoherent 
(phase insensitive) combination of antenna signals achieves the highest event rate per beam, and that 
50-100 MHz processed bandwidth is sufficient for extragalactic searches with SKA Phase 1; the gain 
in event rate from using the full available bandwidth is small. Greater system flexibility will enable 
more effective searches, but need not drive the top-level system requirements beyond those already 
proposed for the SKA. The most appropriate search strategy depends on the observed sky direction 
and the source population; for SKA Phase 1, low frequency aperture arrays tend to be more effective 
for extragalactic searches and dishes more effective for directions of increased scatter broadening, such 
as near the Galactic plane. 
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1 Introduction 

The high time resolution (HTR) universe has been rel- 
atively poorly observed at radio wavelengths and open- 
ing up this parameter space is a key science driver of 
the Square Kilometre Array (SKA), via 'Exploration of 
the Unknown' (Wilkinson et al. 2004). We define fast 
transients as impulsive, singly occurring or intermit- 
tent signals, emitted from high energy density events; a 
search for such events assumes an observed pulse width 
less than the normal correlator averaging time of a few 
seconds. In this context, pulsars can be classed as pe- 
riodic fast transients. Other known transients include 
giant pulses, magnetars and rotating radio transients 
(Cordes 2009; Macquart et al. 2010b). 

Interest in exploring HTR parameter space is grow- 
ing; Tabic 1 lists some searches with existing and fu- 
ture telescopes. The SKA will provide at least an order 
of magnitude improvement in sensitivity over all these 
telescopes and a larger field of view (FoV) than most, 
as shown in Figure 1 of Macquart et al. (2010b). Both 
FoV and sensitivity contribute to the expected rate 
of event detection; employing the analyses in Cordes, 
Lazio, & McLaughlin (2004) and Macquart (2011), this 
yields at least one or two orders of magnitude improve- 
ment in event rate for the SKA. Table 1 compares 
Phase 1 of the SKA (SKAl) to other experiments, us- 
ing the calculations developed in this paper. 

Different search strategies are needed for fast tran- 
sients and pulsars, despite some common HTR require- 
ments. Pulsar surveys proposed for the SKA (e.g. Smits 
et al. 2009; Cordes 2009) take advantage of pulsar pe- 
riodicity to improve sensitivity. They involve a compu- 



tationally expensive, systematic survey of the Galaxy, 
where the observer 'drives' the telescope. The figure of 
merit (FoM) typically used to determine the effective- 
ness of such a survey measures the speed at which an 
area of sky is surveyed to a certain sensitivity (Cordes 
2009). However, this FoM does not consider the num- 
ber of events detectable in a volume of sky, nor does it 
give any weight to the processing cost of sampling the 
sky or searching the data. 

The main goal of at least first-generation fast tran- 
sients searches is to maximise the number of events 
detected in a survey. The expected rate of event de- 
tection depends on, amongst other factors, the search 
strategy employed on the telescope. Given each search 
strategy has a different processing cost, event rate per 
unit cost (T^-t-ost-i) is a more comprehensive FoM than 
survey speed. Put simply, one search strategy may 
have a higher total rate of detection than another, but 
the processing cost of the strategy also needs to be con- 
sidered. This is especially important for SKAl, con- 
sidering transients detection will only be carried out 'if 
it can be done with minimum additional cost or effort' 
(Dewdney et al. 2010). 

Because signal and search processing costs are ar- 
chitecture specific, we use a new FoM, event rate per 
beam formed and searched (7?.boam-i) to generalise the 
problem and parametrise the effectiveness of a search 
strategy. It is based on the rate of transient events 
detectable in a volume of sky as discussed in Mac- 
quart (2011), although it could similarly be applied to 
a FoM for surveying sky area. It assumes T^-t-ost^i cx 
^boam-ij which Is Valid when cost increases linearly 
with the number of beams (independent FoVs) formed 
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Table 1: Radio searches of the high time resolution universe and a comparison of event rate per beam. 



Experiment^ 


Telescope 


^centre 


Av 


Max. 


^bcam 1 


Max. 




and Status 


(MHz) 


(MHz) 


baseline 
(km)*" 


(normalisedj 


beams 
available 


Archival searches'^ 


Parkes 


N/A 


- 


- 


- 


- 


Fly's eye fast radio transient 


ATA 


1420 


210 


N/A 


10-3 (fly) 


42 


search" 


(completed) 












High Time Resolution Universe 


Parkes 


1352 


340 


N/A 


10"2 


13 


Pulsar Survey^ 


(operational) 












Pulsar ALFA (PALFA) 


Arecibo 


1440 


100 


N/A 


10"2 


7 


Survey^ 


(operational) 












V-FASTR^ 


VLBA 


1400 


64 


6000 


10-2 (inc.) 


1 




(operational) 










r 


LOFAR Transients Key 


LOFAR 


120 


32 


< 100 


IQ-^ (inc.) 


Science Project' 


(in progress) 








10-4 (coh.) 
10-2 (inc.) 


thousands^ 


Commensal Real-Time 


ASKAP 


1400 


300 


6 


36 


ASKAP Fast-Transients 


(planned) 








10-s (coh.) 


N/A 


(CRAFT) Surveyj 














Effelsberg Northern Sky Pulsar 


Effelsberg 


N/A 










Survey'*^ 


(planned) 












SKAl AA-low 




260 


380 


200 


1 (inc.) 
10-1 (coh.) 


hundreds^ 
thousands^ 


SKAl low band dishes 




725 


550 


200 


1 (inc.) 
10-2 (coh.) 


1 

thousands^ 



^Only experiments within SKAl frequencies (70 MHz ~ 3 GHz) are listed. Pulsar surveys insensitive to 
single pulses are excluded. N/A is not applicable or information not available. 
''For event localisation using triggered buffer. 

■^Order of magnitude estimation as per Equation 19, normalised to the incoherent combination of SKAl 
low band dishes. For radio telescope arrays, the calculation is for fly's eye (fly), incoherent combination 
(inc.) or coherent combination (coh.), see Section 2. A flat spectrum and no scatter broadening is 
assumed. 

'^Lorimer et al. (2007); Burke-Spolaor k Bailes (2010); Keane et al. (2011). 
"Siemion et al. (2011). 
fKcith ct al. (2010). 

sCordes et al. (2006); Deneva et al. (2009). 
^Wayth et al. (2011). 

'Dutch LOFAR as in Hcssels ct al. (2009); van Lccuwen & Stappcrs (2010). More scenarios are discussed 
in Stappers et al. (2011). 

JMacquart, Hall, & Clarke (2010a); Macquart et al. (2010b). 
'^^https : / / imprs-docs .mpif r-bonn .mpg . de/?p=45. 

Limited by available beamformer processing and data transport. 



and searched. This is true for first-order beamforming 
and data transport costs for the SKA (Chippendale, 
Colegate, & O'SuUivan 2007; Faulkner et al. 2010). 
The search costs also increase linearly because each 
beam signal needs to be searched individually; we con- 
sider the efficiency gain from using a single processing 
unit to process multiple beams to be a second-order 
effect. 

For a given search strategy, T^-bg^jn-i parametrises 
the choice of receptor (antenna), the performance, cost 
and efficiency of the signal combination mode and tran- 
sients search system, and the observed sky direction. 



Because the cost of data storage is prohibitive^ , searches 
are conducted in real-time and data from candidate 
events recorded for subsequent verification and anal- 
ysis. To reduce processing costs, search strategies al- 
ternate to those proposed for pulsar surveys can be 
considered. A pulsar survey need only visit each vol- 
ume of sky once. In contrast, one volume of sky is 
considered as likely as another to contain transients 
(Cordes 2009), and each time a volume is re- visited, 

^For example, ~ 20 GB per station beam per second 
would need to be written to storage for SKAl AA-low 
(380 MHz bandwidth x 2Nyquist X 4 bits X 50 stations). 
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there is new detection potential. This independence 
from sky direction enables lower cost commensal sur- 
veys for fast transients (piggy-back surveys where fast 
transients are not the primary telescope observation). 

Understanding the most effective way to combine 
the signals from the antennas is also important. A ra- 
dio telescope array produces images by correlating the 
antenna signals and averaging the output over a few 
seconds to reduce subsequent processing costs. To ob- 
tain a time resolution of order milliseconds or higher 
for fast transients detection, alternatives to 'fast imag- 
ing' are currently required and various signal combi- 
nation modes must be considered. There are trade- 
offs between a highly sensitive mode with small FoV 
(such as the coherent combination of antenna signals) 
and less sensitive modes which cover more of the sky 
(incoherent combination, subarraying and 'fly's eye'). 
These trade-offs are also influenced by the spatial den- 
sity of antennas in the array. Furthermore, additional 
processing capability enables multiple beams to be formed 
and searched, re-using the array collecting area to some 
extent. The event rate per beam also depends on ob- 
serving frequency and bandwidth. 

The SKA is being designed and constructed in two 
phases, SKAl and SKA2, where the first is a subset 
of the second (Garrett et al. 2010). This paper sets 
out a high-level 'use case' for searching for fast tran- 
sients with SKAl receptors: low frequency aperture 
arrays and low band dishes (Section 2). It outlines 
a basis for comparing survey strategies (Section 3) 
and undertakes a detailed analysis of the effects of 
receptor choice, signal combination modes, sky direc- 
tion, observing frequency and bandwidth on the event 
rate per beam (Sections 4 and 5). These effects are 
summarised (Section 6) and specific recommendations 
for SKAl are made (Section 7). This paper extends 
the work done by the International Centre for Ra- 
dio Astronomy Research (ICRAR) for the Commensal 
Real-Time ASKAP Fast-Transients (CRAFT) survey 
(Macquart et al. 2010a,b). CRAFT is one of the sur- 
vey science projects planned for the Australian SKA 
Pathfinder (ASKAP), a designated SKA precursor in- 
strument. The analysis is applied to the SKAl system 
description to determine optimal search strategies, but 
the method is equally applicable to SKA2 and other 
radio telescope arrays. 

2 Fast transients search use case 

For context, this section outlines a generic fast tran- 
sients search use case. Because the data volumes are 
too large to store cost-effectively, the search for tran- 
sient events is conducted in real-time on a data stream 
which is a continuous observation of the sky. How- 
ever, a rolling buffer recording a small period of the 
data as it is observed allows candidate events contain- 
ing potential fast transients detections to be saved and 
further processed off-line. The specific implementa- 
tion of a fast transients search depends on the tar- 
get or expected source population and the cost and 
performance factors of components in the processing 
pipeline; see Macquart et al. (2010b), Stappers et al. 



(2011) and Wayth et al. (2011). Figure 1 shows a 
generic fast transients pipeline and the signal process- 
ing steps (actions) in this pipeline are as follows: 

Signal reception: Radio signals are collected by 
SKAl receptors: low frequency (70-450 MHz) aper- 
ture arrays (AA-low) and a dish array equipped with 
low (0.45-1 CHz) and high (1-2 CHz) band single pixel 
feeds. Only the low band feed on the dishes (hereafter 
referred to as low band dishes) is discussed in this pa- 
per, because it generally achieves a higher event rate 
than the high band feed, given the assumptions in Sec- 
tion 4.1. 

Station beamforming: The complex signals of 
many AA-low elemental antennas are combined vec- 
torially to form one or more station beams, as shown 
in Figure 2. The FoV of each station beam is deter- 
mined by the station diameter and is similar to a single 
dish of the same diameter. The output station beam 
can then be processed in the same way as a dish sig- 
nal. Though not considered here, signals from a group 
of dishes (instead of elemental antennas) can similarly 
be coherently combined into dish station beams. Note 
that both dish and AA stations are simply subarrays 
of coherently combined antennas. 

Signal combination: The signals from the dish 
or station beams can be combined coherently, incoher- 
ently or not at all. These signals are then searched for 
fast transients. Figure 2 shows the signal combination 
modes and resultant beam patterns; they are further 
discussed in Appendix A. Incoherent (phase insensi- 
tive) combination sums the detected signals (powers) 
from antennas pointing in the same direction. Coher- 
ent combination of antennas forms a phased or tied 
array beam — voltages measured at each antenna are 
aligned in phase towards a specific direction on the sky, 
in a manner similar to station beamforming. Smaller 
groups of antennas — subarrays — can be incoherently 
combined and each subarray pointed in a different di- 
rection. The extreme of this is the so-called fly's eye, 
where every antenna is pointed in a different direction. 

Event localisation and the spatial discrimination of 
astronomical signals from radio frequency interference 
(RFI) is possible for coherent combination and, us- 
ing buffered voltages, for incoherent combination and 
subarrays of three or more incoherently combined an- 
tennas. Multiple beams (incoherently or coherently 
combined) can be used to discriminate RFI, where a 
candidate event in most or all beams indicates the pres- 
ence of RFI. 

Dedispersion processing: The signals pass through 
a cosmic medium of unknown dispersion measure (DM). 
This means that the detection needs to be trialled for 
many DMs, each of which has a computational cost. 
The DM range to be trialled depends on the location 
on the sky. Clarke et al. (2011) discusses dedispersion 
for SKAl in detail. 
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Figure 1: High-level flow diagram for a generic fast transients pipeline, for SKAl receptors. Rounded 
boxes are signal processing actions, rectangles describe the information flow. The solid vertical line is the 
data spigot point in the signal chain (see text for details). 



Event detection: An event detection algorithm 
needs to be applied to the signal from each trial DM, 
where optimal detection is achieved with an appropri- 
ate matched filter (Cordes & McLaughlin 2003). 

Store in rolling buflfer: The digitised voltages 
from the dishes or stations are stored in a circular 
memory (rolling) buffer. In the case of a candidate 
event, the data from the buffer can be saved to an- 
other location (dumped) and processed off-line. The 
amount of memory required in the buffer depends on 
the sampling rate, sample size and the expected maxi- 
mum (dispersed) pulse duration. The maximum pulse 
duration is a function of the range of frequencies to be 
captured and varies linearly with the maximum DM 
to be trialled. For a maximum DM of 3000 pc cm"'' 
and a bandwidth of a few hundred MHz, a buffer of 
order tens of seconds is required for dish frequencies 
and possibly tens of minutes for lower frequencies. 

Buflfer dump and ofF-line processing: On re- 
ceipt of a trigger, the buffer will dump the original 
voltage data to storage for off-line processing, which 
could include RFI filtering, analysis of the candidate 
detection and correlation of the dish or station beams 
for source localisation and imaging. 

Commensal and targeted surveys: A com- 
mensal survey greatly increases observation time by 
conducting the survey in parallel with normal telescope 
operations. It is passive; it uses dish or station beam 
signals from the primary user observation, placing lit- 
tle extra demand on the telescope. Such a survey is 
suitable for extragalactic searches, given the informa- 
tion about the population of such fast transients is 
not fcnown a priori; hence one direction on the sky 
is as good as another. To observe specific areas of 
the sky, such as the Galactic plane and nearby galax- 
ies, a targeted transients survey (which is the primary 



user observation) may be required (e.g. van Leeuwen 
& Stappers 20f0). 



Data spigot A data spigot to the dish and sta- 
tion beam signals is useful for transients surveys, es- 
pecially those which are commensal. If the signal chain 
is considered to be the signal path from the antennas 
of a radio telescope array to the correlator, a spigot 
defines a point in the signal chain where users can 
tap off data via a well defined interface. The spigot 
for fast transient searches may output either coher- 
ent (phase-preserved) data at high rates or, alterna- 
tively, incoherent data where the dish or station beam 
voltages are squared and integrated to a time resolu- 
tion of order milliseconds to reduce the data rate and 
subsequent dedispersion processing load. The latter 
approach is being taken by CRAFT to access beams 
from the ASKAP beamformer (Macquart et al. 20f0a). 
Similarly, searching the integrated signals from the 
dish or station beams which have been incoherently 
combined is a low cost option for commensal surveys 
with SKAl. 

The solid vertical line in Figure 1 shows the point 
in the fiow diagram where the spigot for fast tran- 
sients would need to exist to enable the signal com- 
bination modes in this paper. The pipeline after the 
spigot point is not part of the normal imaging mode 
of the telescope; the post-spigot pipeline may be im- 
plemented internally or with user-provided processing. 
An example of processing being implemented inter- 
nally is the 'non-imaging processing' for pulsar ob- 
servations with SKAl (Dewdney et al. 2010); the ap- 
proach being taken by CRAFT is an example of user- 
provided processing. Note that it is conceivable that a 
spigot to the AA-low elemental antenna signals could 
also exist, but the data rates make this option pro- 
hibitively expensive for SKAl. 
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Single element: Incoherent Coherent Incoherently 

AA station or dish combination combination combined subarrays 




Figure 2: Signal combination modes, resultant beam patterns, and beam terminology for dishes and 
aperture array (AA) stations (beam sizes not to scale). 



3 Survey strategy: Maximis- 
ing survey speed and min- 
imising cost 

Figures of merit quantify the effect of ahering the vari- 
able parameters of a problem. A simple FoM to mea- 
sure the cost effectiveness of a fast transients search 
strategy is the detected event rate per beam searched 
('^beam-i)i ^ proxy for cost in the absence of suffi- 
ciently accurate design and cost information. We want 
to optimise for a high value of TZy^^^^-i , although the 
total event rate for all beams (TZ,^) must also be high 
enough to be of scientific benefit and open new vol- 
umes of parameter space. The qualitative advantages 
of a strategy, which cannot be captured in a FoM, must 
also be considered. 

The detected event rate is effectively a survey speed. 
Smits et al. (2009) present a frequency-dependent FoM 
for survey speed (SSFoM) for dishes. It is based on sur- 
veying an area of sky, thus SSFoM is linearly propor- 
tional to FoV and sensitivity squared. In this paper, 
the equivalent SSFoM is the rate of transient events de- 
tectable in a volume of sky and is linearly proportional 
to FoV and sensitivity to the power of 3/2. It draws 
on event rate calculations from Macquart (2011); the 
derivation is shown in Appendix B. 

The event rate is given by 

7?.,. = p^^^^Fmax events s"\ (1) 
47r 

where Vmax is the maximum volume out to which an 
object is detectable and Q,proc is the processed FoV, 
which is the product of the number of beams formed 
(A''boam) and the FoV of each beam. Thus the event 
rate per beam is 

^bcam-l = 1^ • (2) 

*bcam 

A^beam may be the number of station (Ab-o) or array 



(Ab-arr) beams. For subarrays, it is the product of 
the number of subarrays and station beams formed in 
each subarray (AsaAi,_o) • 

We describe an extragalactic survey as a search for 
a homogeneously distributed population of isotropi- 
cally emitting fast transients of fixed intrinsic lumi- 
nosity. For such a population, the event rate is 

= ^Pf^p- (f ) ' (4J;-) ' eventss-\ (3) 

where p (events s~^ pc~^) is the event rate density, 
d (Jy pc^) is the intrinsic luminosity of the population, 
Wi is the intrinsic pulse width, W is the observed pulse 
width and Smin is the minimum detectable flux density 
of the telescope for an integration time oir = Wi. The 
(Wi/VK)^'''' term approximates the loss in signal-to- 
noise (S/N) due to pulse broadening (see Appendix B). 
The frequency dependence of TZv is discussed in more 
detail in Section 4.3. 

Macquart (2011) shows that the event rate is pro- 
3/2 

portional to f2proc and 5^^^;^^ for an extragalactic pop- 
ulation with a luminosity distribution which follows 
a power-law or lognormal distribution. Although the 
actual event rate depends on the luminosity distribu- 
tion, the proportionality still holds for Slproc and 5min, 
which is sufficient to compare telescopes and their sig- 
nal combination modes. For fast transients searches 
within the Galaxy, scatter broadening due to multi- 
path propagation in the interstellar medium makes the 
event rate per beam dependent on frequency and di- 
rection. Although comprehensive direction-dependent 
modelling is beyond the scope of this paper, the loss 
in sensitivity due to scattering is incorporated in the 
W term in Equation 3 and modelled for some repre- 
sentative sky directions. This first-order analysis gives 
an indication of how SKAl will perform as a function 
of frequency; see Section 5.3 for further discussion of 
Galactic objects. 
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4 Modelling event rates 

Event rate per beam is a simple metric to model fast 
transients event rates in the cost constrained environ- 
ment of SKAl. This section shows how signal com- 
bination mode and filling factor (antenna spatial den- 
sity) affect Ti-iiciiin-^ fo'" * radio telescope array. The 
dependence on frequency and the effectiveness of search- 
ing large bandwidths are also considered. The mod- 
elling is specifically applied to the SKAf receptors, for 
source populations whose intrinsic luminosity either 
does not vary with frequency or varies with 
value typical of the pulsar population (Lorimer et al. 
1995). This paper uses the system description of SKAl 
(Dewdney et al. 2010) to make trade-offs — the relevant 
details are given in Table 2. However, the complete 
system design for SKAl is still under development and 
subject to a decision-making process involving trade- 
offs and performance and cost optimisation. 

4.1 Assumptions 

The trade-offs in this section make the following sim- 
plifying assumptions: 

• The population of fast transients is homogeneously 
spatially distributed and of fixed intrinsic lumi- 
nosity. 

• A matched filter is used to detect the dedis- 
persed, but scatter broadened pulse (as per Cordes 
& McLaughlin 2003). 

• The effects of scintillation on source intermit- 
tency and optimum search bandwidth are ig- 
nored. 

• The dedispersion processing system does not con- 
tribute to pulse broadening (see Appendix C 
for a description of these instrumental contri- 
butions). 

• The intrinsic pulse width is 1 ms. Shorter dura- 
tion pulses would be more sensitive to S/N loss 
due to pulse broadening, longer duration pulses 
less sensitive. 

• Events are broad-band such that the intrinsic 
spectral bandwidth of the pulse is greater than 
the processed bandwidth. Thus all channels across 
the band contain contributing signal. 

• The beam has constant (maximum) sensitivity 
between the half-power beamwidth points, and 
zero sensitivity outside of that. 

• Beamformer calibration costs are not considered. 

• A time to frequency domain transformation (chan- 
nelisation) and cross-correlation 'FX' correlator 
is used, as opposed to other correlator topolo- 
gies such as 'XF'. It is the most cost-effective 
architecture for the SKA, and allows other sig- 
nal processing actions, such as beamforming and 
RFI excision, to be done efficiently (Hall et al. 
2008). 

• The processing cost of forming and searching a 
beam is independent of frequency, bandwidth 
and signal combination modes. In practice, lower 



Table 2: SKAl system details. 

Low frequency aperture arrays (AA-low) 
Aperture 

Frequency range^ 70-450 MHz 

Station diameter (Dq) 180 m 

Number of stations (iVo) 50 

Number of antennas (A^a) H 200 per station 

Station beam taper (JCq) 1.3 

Dense-sparse transition 115 MHz (2.6 m) 

(^transition) 

Array configuration regions'' 

Core (radius<0.5 km) ^ 50% (25 stations) 

Inner (l<radius<2.5 km) 20% (10 stations) 

Mid (2.5<radius<100 km) ~ 30% (15 stations) 
Core filling factor 0.81 
Performance 

Trcvr 150 K 

Bandwidth per beam (Aj/) 380 MHz 
Single pixel feed dishes 

Aperture 

SKA2 dish frequency 0.3-10 GHz 
capability 

Parabolic dish diameter 15 m 
(^o) 

Number of dishes (^o) 250 

Total physical aperture 44 179 m^ 

Dish illumination factor 1.15 
(/Co) 

Array configuration regions'' 

Core (radius<0.5 km) - 50% (125 ant.) 

Inner (0.5<radius<2.5 km) - 20% (50 ant.) 

Mid (2.5<radius<100 km) - 30% (75 ant.) 

Core filling factor 0.03 
Antenna RF system'^ 

Feed/LNA low band 0.45-1.0 GHz 

Bandwidth (Az^iow) 0.55 GHz 

Feed/LNA high band 1.0-2.0 GHz 

Bandwidth (Ai^high) 1-0 GHz 
Performance 

Antenna/feed efficiency'^ 70% 

Average Tgys in low band" ^40 K 

Average Tgys in high band ^30 K 

^ Single dual polarization antenna over frequency 
range. 

^ Fractional number in each region. 

One dual polarization feed available at a time. 

Average over frequency. 
° Higher at the low frequency end of this band. 
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frequencies (where the maximum dispersed pulse 
duration is longer) and larger bandwidths will 
increase processing costs; the magnitude of the 
increase is specific to the processing architecture 
and the effect of sky direction on the DM range 
to be trialled. 

4.2 Signal combination mode com- 
parisons 

This section compares the per beam event rate for in- 
coherent and coherent combination and fly's eye (Fig- 
ure 2 and Appendix A) and applies these results to 
SKAl AA-low and low band dishes. We refer to a sin- 
gle dish or A A station as an element, designated with 
the subscript 0. The combination of an array of No 
elements may refer to the number of elements in the 
total array, or some subset of the total array (e.g. the 
SKAl core region in Table 2). 

Incoherent combination of an array of A'o elements 
increases the sensitivity by a factor of \/TVo over a sin- 
gle element while retaining its FoV, ilo- Forming A^b-o 
station beams linearly increases the FoV. To even fur- 
ther increase the FoV, A'^sa subarrays can be incoher- 
ently combined. Each subarray is pointed in a differ- 
ent direction, increasing the FoV by a factor of Aga but 
only increasing the sensitivity of the array by a factor 
of \/Ao/sa over a single element, where Ap/sa is the 
number of elements per subarray. Fly's eye pertains 



to the case N, 



O/sa 



= 1. 



Sensitivity of the coherent combination of an array 
of Ao elements is higher than incoherent combination 
and subarraying; it increases proportional to A^o- How- 
ever the FoV of the array beam, fiarr, is much smaller; 
it is proportional to D~^^, where Darr is the diameter 
of the array of elements being combined. The FoV can 
be linearly increased by forming Ab-arr array beams. 

Applying these relationships to Equation 3 gives 
the total event rate for each signal combination mode: 



^1 / WiApoiAt/r \i / CiA^o ' - 



M, 



an event rate FoM equivalent to one incoherently com- 
bined beam for ASKAP. We modify this analysis and 
apply it to TZv. 

Following Cordes (2009), we define the number 
of pixels (Apix) as the maximum number of indepen- 
dently pointed, coherently combined array beams which 
can be formed within the FoV of a single element beam. 
It is frequency independent, and given by 



'Ab-ot^oAo'''* 



M = <N, 



b — arr^arr Aq 



3/2 



Incoherent combination 
Coherent combination 



AsVAb-of^oAQ'''* Subarraying 



(4) 



where Apd is the number of polarisations summed, Ai/ 
is the processed bandwidth, r is the post-detection in- 
tegration time (which also defines the time resolution 
of the observation) , A^o is the effective area of an ele- 
ment (dish or station), a is the S/N ratio required for 
event detection and Tgys is the system temperature. 



4.2.1 Filling factor efficiency 

The coherent combination mode is more effective if 
the dishes or stations are closely spaced, thus having a 
higher filling factor. In this case, the same number of 
elements are being combined, but the array beam FoV 
is larger. D'Addario (2010) considers the number of 
coherently combined array beams required to achieve 



no 

f^ar 



ICo -Dar 
^arr Do 



(5) 



where ICo and ICmi are the element and array beam 
tapers respectively. 

A measure of the effectiveness of the coherent com- 
bination mode is the number of array beams which 
need to be formed and searched to achieve a coher- 
ent combination event rate (7?.coh) equal to incoherent 
combination (7?.inc). From Equation 4, 7?.coh = 7?.inc 
gives 



Ah 



Ab-o»oAg/^ 

na„(7?Ao)3/2 

Ab^oApix 



(6) 



3/4 



where rj is the fraction elements in the array which are 
coherently combined, out of a total No- For example, 
77 = 0.5 if only the elements in the SKA core are co- 
herently combined while those in the total array are 
incoherently combined. 

Achieving the highest possible event rate is desir- 
able, but this must be tempered by the cost of search- 
ing multiple beams. The relative event rate per beam 
depends on the array filling factor and is simply the 
inverse of Equation 6 when Ab-o ~ 1: 



nc beam ^ • 



(7) 



For incoherently combined subarrays (when all Aq el- 
ements are formed into subarrays). 



7? 1 — N^'^^'^V 

'^sabcam ^ -'^sa '^inc bi 



cam ^ ' 



(8) 



where Asa subarray beams are searched. 



4.2.2 Coherent combination for a fully filled 
array 

Assuming a best case scenario of an array entirely filled 
with stations of equal diameter and tCo = A^arr, 



Ap 



Ao. 



(9) 



Substituting into Equation 6 and for rj = 1, the the- 
oretical minimum number of array beams required so 
that TZcoh = 7?.inc is 

Ab-arr = Ab-oAt}/* (10) 

and the relative event rate per beam (Equation 7) is 



^cohl 



1 —^'0 'Mncbcam-l- 



(11) 
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This estimation is optimistic towards coherent combi- 
nation, given that it is not physically possible to en- 
tirely fill a circular array with circular stations. Re- 
gardless of this, for A'^o > 1, the incoherent combina- 
tion will always achieve a higher event rate per beam 
searched than coherent comhinatton, and this differ- 
ence increases with No. These results are frequency 
independent. 

4.2.3 Signal combination modes for SKAl 

A flexible processing system allows various signal com- 
bination modes; Table 3 compares the event rate for 
incoherent and coherent combination and fly's eye, for 
SKAf receptors. The number of array beams required 
such that 7?.coh = Tiinc is calculated for coherent com- 
bination. The event rate per beam formed and searched, 
relative to T^-in^-bcam-i i is calculated for all modes. Three 
coherent combination modes are shown: elements in 
the core region, inner and core region and the total ar- 
ray. Of these, using the elements in the core achieves 
the highest event rate per beam, due to the higher den- 
sity of collecting area. For this reason, the core will be 
the only coherent combination mode further analysed 
in this paper. 

The coherent combination of the AA-low core re- 
quires approximately five array beams to equal the 
event rate of a single beam of the incoherently com- 
bined total array. Although the fiy's eye mode achieves 
a higher event rate by a factor of SO^'^* = 2.6, the rela- 
tive rate per beam is less than the coherently combined 
core. If the dense packing is not achievable (Dewdney 
et al. 2010), the number of array beams required will 
be higher. Indeed, using optimisations from Graham 
et al. (1998), the optimal packing of 25 congruent cir- 
cles of diameter 180 m in a circle results in a minimum 
core diameter of 1036 m. This is larger than the 1000 
m diameter in Dewdney et al. (2010) and excludes any 
spacing that may be required for infrastructure. 

Due to the lower filling factor, 157 beams formed 
from the coherent combination of the dishes in the core 
are required to equal the event rate of the incoherently 
combined array. The event rate of the fly's eye mode 
is higher by a factor of 250^^* = 4 and its relative 
rate per beam is higher than the coherently combined 
core. Although not shown, values for subarrays lie 
between fiy's eye and incoherent combination modes, 
and depend on the number of elements per subarray. 

4.3 Frequency dependence 

The event rate has a frequency dependence on luminos- 
ity, minimum detectable flux density, FoV and scatter 
broadening, designated by subscript i/: 



Wi 



\ 47r5min„ 



(12) 



It is important to understand this frequency depen- 
dence because of the wide fractional bandwidths of the 
SKA. Looking at each of these dependencies in turn: 

• The processed FoV depends on the number of 
beams formed and whether they are formed in- 



coherently or coherently, but either way is pro- 
portional to v~^: 



cK. 



(13) 



where c is the speed of light. For incoherent 
combination, A'^beam is the number of station 
beams formed, K. is the feed illumination factor 
or station beam taper and D is the diameter of 
the dish or station. For coherent combination, 
Nheaia is the uumbcr of array beams formed, K, 
is the array beam taper and D is the diameter 
of the array (see Appendix A). 

• Pulses are broadened due to scattering. The 
broadening time Td depends on the path through 
the Galajcy to the observer, and scales as Td oc 
^-4.4 ^QQj.(jgg ^ Lazio 2002). The observed pulse 



duration is given by 



(14) 



See Appendix C for further details. 

• Because we are looking for an unknown popula- 
tion, we do not know how the luminosity varies 
with frequency. However, spectral indices for 
pulsars have been measured. We use ^ = —1.6 
(Lorimer et al. 1995), a value typical of the pul- 
sar population, such that 



(15) 



where Co is luminosity at reference frequency i/q . 
For comparison, we also consider a fiat spectrum 
population (^ = 0). 

• Smin„ is a function of Tsys^ and Aco^ : 



(16) 



For aperture arrays, the effective area of a sta- 
tion is approximately 



^ ^transition 



-^a X V > /^transition- 



(17) 



The AA-low system temperature is the sum of 
the receiver noise and an approximation to the 
sky temperature: 



Tsys TycvT ~\~ 60 



(18) 



A graphical breakdown of the frequency dependencies 
of AA-low is shown in Appendix D. 

The dependence of event rate on frequency for three 
representative sky directions and spectral indices of 
^ = or —1.6 is shown in Figure 3. The normalised 
^boam-i is plotted at 1 MHz intervals for centre fre- 
quency v and processed bandwidth Ai^ = 1 MHz, span- 
ning the SKAl system description frequency range of 
AA-low and low band dishes (70 MHz < v < 1000 
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Table 3: Relative event rates for SKAl receptors and select signal combination modes 



Signal combination mode 


Input parameters'^ 




Calculated values 






Darr (km) 


No 


A^pix 


-^b-arr 

required 


Relative TZ 
bcam"^ 


per 


AA-low receptors 












Incoherent: total array 




50 






1 




Coherent: core 


1 


25 


31 


5A'b-o 


2.15 X 10 


-1 


Coherent: inner + core 


5 


35 


77z 


70A^b-o 


1.43 X 10 


-2 


Coherent: total array 


200 


50 


1.23 X 10^ 


6.57 X lO'^iVb-o 


1.52 X 10" 


-5 


Fly's eye: total array 




50 

subarrays 






5.32 X 10" 


-2 


Low band dish receptors 












Incoherent: total array 




250 






1 




Coherent: core 


1 


125 


3.49 X 103 


157 


6.37 X 10- 


-3 


Coherent: inner + core 


5 


175 


8.72 X 10^ 


2.37 X 10^ 


4.22 X 10" 


-4 


Coherent: total array 


200 


250 


1.40 X 10^ 


2.22 X 10^ 


4.50 X 10" 


-7 


Fly's eye: total array 




250 

subarrays 






1.59 X 10" 


-2 



From Dcwdney et al. (2010). 
^ Number of array beams required for an event rate equivalent to incoherent combination. For dishes, 
iVb-o = 1. 

Relative to the incoherent combination event rate per beam. 



MHz). The relative event rate per beam between sig- 
nal combination modes (the rightmost column of Ta- 
ble 3) still applies as a multiplicative factor to the data 
in Figure 3, regardless of sky direction and frequency. 

The simplest case to consider is a search for extra- 
galactic fast transients (solid line in Figure 3). In this 
case there is no sensitivity loss due to scatter broaden- 
ing; VF S> Td is assumed. Cordes & McLaughlin (2003) 
find that for a given sky direction, the scatter broad- 
ening of an extragalactic source will be approximately 
six times the broadening from the Galaxy alone. This 
assumes equal scatter broadening in the host galaxy, if 
there is one, and no contribution from the intergalactic 
medium. For directions away from the Galactic plane 
where Td is low and for an intrinsic pulse width of 1 
ms, the exclusion of scatter broadening is a reasonable 
first-order assumption. 

A first-order analysis of a search for Galactic tran- 
sients is possible by invoking the simplifying assump- 
tions listed in Section 4.1. We calculate relative event 
rates taking into account estimates of scatter broaden- 
ing at a distance of 30 kpc and frequency of 1 GHz from 
the NE2001 model of Cordes & Lazio (2002), where the 
broadening scales as rd oc A distance of 30 kpc 

determines the maximum broadening due to interstel- 
lar scattering for that direction; broadening is less at 
shorter distances. Two representative sky directions 
for Galactic transients are: 

• above the Galactic plane: Td = 0.03 ms {DM = 

215pccm"^) at 1 = 300, b = 10 

• on the Galactic plane: Td = 1.78 ms {DM — 
628pccm"^) at 1 = 300, b = 0. 



For these directions, the normalised 7?.boam-i in Fig- 
ure 3 shows how increased scatter broadening reduces 
the event rate at lower frequencies. In this first-order 
analysis, the effect of scatter broadening on event rate 
is independent of the signal combination mode chosen; 
the frequency-dependent (VFi/IF)'^''^ term in Equation 4 
is common to all modes. The factors to consider for a 
more detailed analysis of Galactic transients searches 
is outlined in Section 5.3. 



From Figure 3, the incoherent combination of dish 
signals would be most efficient for searching for ex- 
tragalactic sources with a low spectral index. AA-low 
would be more efficient for such sources with high spec- 
tral indices. For a Galactic population, the preferred 
receptor depends on the amount of scatter broadening. 
For directions near or on the Galactic plane, low band 
dishes show a higher Ttb^a.m-^ than AA-low. However, 
the difference between the two receptors is less than 
an order of magnitude for steep spectrum sources near 
the Galactic plane. For steep spectrum sources, even 
those on the Galactic plane, the higher event rate is at 
the lower end of the frequency band of each receptor. 
For shallow spectrum sources near or on the Galactic 
plane, there is no strong maximum within a receptor 
frequency band. For dishes, the maximum may be at 
the low (450 MHz) or high (1 GHz) end of the fre- 
quency band, for small or large Td respectively. For 
AA-low, the event rate quickly reduces for frequencies 

below t/transition = 115 MHz. 




Figure 3: Normalised event rate per beam for the incoherent combination of the total array for Ai^ = 1 MHz 
and a source spectral index of -1.6 (left) and (right), at three representative sky directions. The rate 
for the coherent combination of the elements in the array core is less by a factor of approximately five for 
AA-low and 157 for the low band dishes. The post-detection integration time equals the intrinsic pulse 
width: T = Wi = 1 ms. Scatter broadening and DM is calculated ai — 1 GHz from NE2001 (Cordcs & 
Lazio 2002) for a distance of 30 kpc and scaled using Td oc v~^-^. Data is normalised to 'R-hca.m-^ = 1 at 
vq for the extragalactic case. The aperture array dense-sparse transition frequency is at 115 MHz. 



Table 4: Normalised extragalactic event rate per beam'' for the full SKAl receptor bandwidth. 



Receptor 


Incoherently combined: 


total array 


Coherently combined: core 




C=-1.6 


C=o 


e=-1.6 e=0 


Low band dish 


989 


262 


6.30 1.67 


{Au = 550 MHz) 








AA-low (Azy = 380 MHz) 


3.13 X 10" 


164 


6.74 X 10^ 35.3 


AA-low and low band dish 


3.15 X 10" 


361 


6.74 X 10^ 35.7 



^ Normalised to 7?.boam-i = 1 = 1 MHz) for the incoherent combination of array elements at vq = 
1 GHz. The actual event rate per beam is 6.94x pC^^x events s^^, where x is the normalised 7?.bcam-i i 
p has units of events s^^ pc^'^ and £o has units of Jypc^. 



4.4 Large processed bandwidths 



iVchAi^ch using 



Processed bandwidth (A;/) is the bandwidth of the as- 
tronomical signal at the fast transients detection sys- 
tem. In radio astronomy, an increase in processed 
bandwidth is usually assumed to produce a ^/~Ku in- 
crease in signal-to-noise. From Equation 4, this would 
produce a Au^^'^ increase in event rate. However, un- 
less the event rate is approximately constant across 
frequency, this increase does not hold for large pro- 
cessed bandwidths. 

To calculate the event rate over a large processed 
bandwidth, the frequency-dependent channel contri- 
butions shown in Figure 3 are summed. Where the 
pulse broadening (due to propagation and instrumen- 
tal effects) is not significant, the event rate is calcu- 
lated numerically for a processed bandwidth of A.u = 



TZav 




(19) 

where A^ch is the number of frequency channels of width 
A//ch- The derivation is shown in Appendix E. Table 4 
makes these calculations for different modes and spec- 
tral indices, for the full bandwidth of AA-low and low 
band dishes. The relative event rate per beam between 
signal combination modes of a given receptor (Equa- 
tions 7 and 8) still hold for large processed bandwidths. 
Strikingly, for the combined AA-low and low band dish 
event rate, the contribution from the dishes is only sig- 
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■ — c = o 






C = -l-6 — 








■ J Dishes (ch 1 


= 450 MHz) 










$ AA-low (ch 1 


= 115 MHz) 









)0l 1 1 1 1 ^ 

100 200 300 400 500 



Number of 1 MHz channels 

Figure 4: Extragalactic event rate, calculated for 
1 to 335 AA-low channels of width 1 MHz, where 
channel 1 is at 115 MHz, and 1 to 550 low band 
dish channels, where channel 1 is at 450 MHz. 
Thick lines show the calculated TZal/ versus pro- 
cessed bandwidth Aiy = Nch MHz, for spectral in- 
dices of and -1.6. The thin line is the expected 
/\j/3/4 ijicrease in event rate over the rate using 
only one channel (7?.ai/=i mhz)- The rate for each 
curve is normalised such that 7?.Aiy=iMHz — 1; the 
slope of each curve applies equally to the incoher- 
ent or coherent combination of the receptor. 



11 



For example, say the threshold is set to 0.5%. Then for 
the ideal case, A'^ch = 150 channels contribute to TZav 
Adding a 151*^* channel will contribute less than 0.5 % 
to 7?.Ai^. For the two cases plotted. Table 5 shows the 
maximum channel for which the improvement in event 
rate over the rate without that channel is greater than 
0.5%. 

Table 5: Maximum number of channels of band- 
width Avch = 1 MHz contributing more than 0.5% 
of the cumulative event rate. For the ideal case, 
the maximum channel number is 150. 



Case Channel number 

^ = -1.6 C = 

AA-low, 47 69 

chl = 115 MHz 

Low band dishes, 87 112 

ch 1 = 450 MHz 



To interpret this table, compare the maximum con- 
tributing channel number when ^ = —1.6. For AA- 
low, 47 channels {Av — 47 MHz) contribute above the 
threshold. For dishes, this is achieved with Av — 
87 MHz, implying that the AA-low processed band- 
width becomes less useful more quickly. This is ex- 
pected, given the steeper spectral dependence of AA- 
low over dishes, shown in Figure 3. 



nificant for the incoherently combined total array case, 
when ^ = 0. 

Given Av = NchAi/ch, plotting 7?.ai/ as a function 
of the number of contributing channels shows the de- 
creasing contribution of higher frequency channels to 
the event rate. This decreasing contribution results in 
the event rate increasing by less than Au^'''^. The only 
exception is in the dense AA regime (< 115 MHz) for 
low spectral indices (—1 < ^ < 0), where the decreas- 
ing contribution comes from lower frequency channels. 

Consider two cases for extragalactic searches, plot- 
ted in Figure 4: the AA-low band from 115 to 450 MHz 
and the low band dishes where the whole 550 MHz 
bandwidth is available. The ideal case of a A//^/* in- 
crease over the Av = 1 MHz event rate is also shown. 
As expected, this plot shows that the maximum band- 
width achieves the highest event rate. However, the 
event rate curve flattens out well before the maximum 
bandwidth, especially for AA-low and also steeper spec- 
trum sources. 

Assuming a limited amount of signal processing is 
available, the following question arises: at what point 
can the processing be more effectively used elsewhere 
in the fast transients pipeline, and how is this quan- 
tified? One method is to arbitrarily set a threshold 
beyond which additional channels contribute very lit- 
tle to the event rate. Channels of increasing frequency 
are included while the following is true: 

^'t;"^'^"'" > threshold. (20) 



5 Discussion 

5.1 Effectiveness of combination modes 

Comparing signal combination modes and array con- 
figuration using event rate per beam (T^-bgam-i) gives 
a frequency independent analysis of the trade-offs for a 
given receptor. Incoherent combination always achieves 
a higher TZy^^^^-i than coherent combination, as shown 
in Section 4.2.2. The difference increases with the 
number of elements or a reduced filling factor. The 
per beam analysis is important because it captures 
the first-order processing costs for each beam. For this 
reason, although fiy's eye gives the highest total event 
rate, it comes at the expense of searching many more 
beams. 

However, the per beam analysis does have some 
limitations. The total event rate TZi, of a search strat- 
egy is 

7?.^ = 7?.bcam-l Abeam, (21) 

where Abeam is the number of beams formed and searched. 
Depending on the signal combination mode used. Abeam 
may be the number of station (A'^b-o) or array (Ab-arr) 
beams, or NaaN-t-o for A^'sa incoherently combined sub- 
arrays. Electromagnetic and signal processing system 
design and cost considerations will put maxima on each 
of these. For example, the number of station beams 
which can be formed within the antenna FoV is lim- 
ited by performance degradation, the diameter of the 
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station and the beamforming processing power avail- 
able. Once the limit on forming more station beams 
is reached, TZi, will eventually become higher for co- 
herent rather than incoherent combination. This is 
because iVpix array beams can be formed within each 
of the A'^b-o station beams (Equation 5), allowing a 
maximum of A^pixA'^b-o array beams to be formed. 

Dishes with single pixel feeds only have a single 
beam (i.e. A'^b-o = 1), so TZi, of the incoherent combi- 
nation of dishes cannot be increased by forming more 
beams. However, if phased array feeds (PAFs) are 
available on dishes, then A'b-o dish beams can be formed, 
increasing the total event rate. 

So we must consider: do we want to form and 
search more array beams, subarrays or incoherently 
combined station beams? This will be influenced by 
station beamforming (for AAs) and array beamform- 
ing costs, search costs and the array filling factor. Al- 
though exact costs are unavailable, we can generalise 
preferred combination as follows: 

• Incoherent combination 

Achieves the highest event rate per beam, mak- 
ing it preferable in most cases. If the cost of 
searching a beam signal is high, incoherent com- 
bination presents a further advantage over other 
modes. 

• Coherent combination 

Requires an array with a high filling factor and 
low beam search costs; it can achieve the high- 
est total event rate if many array beams are 
be formed and searched. For AA-low, the cost 
of forming multiple coherently combined array 
beams must be lower than the cost of form- 
ing and incoherently combining multiple station 
beams. 

• Incoherently combined subarrays and fly's 
eye 

This mode is only preferable when the beam 
search cost is low. A fly's eye mode excludes the 
buffering and source localisation advantages of 
an array, and commensality with most observa- 
tions (see Section 2). Sets of three-element sub- 
arrays counter this problem and could employ 
antennas unused by the primary user observa- 
tion. For example, splitting 24 AA-low stations 
outside the core into 8 of these three-element 
subarrays results in a total event rate approxi- 
mately equal to the incoherent combination of 
50 stations. 

The SKAl AA-low results show a weak preference 
for the incoherent combination of station beams: five 
array beams formed from the coherent combination of 
the AA-low stations in the core are required to equal 
the event rate of a single beam of the incoherently 
combined total array, assuming the stations are very 
closely packed. If multiple station beams are formed 
'at no cost' for normal (imaging) array observing, or 
the cost of array beamforming is high, the advantage 
of incoherent combination is increased. In practice, 
other effects such as RFI mitigation and station and 
array beam quality also need to be taken into account. 



For SKAl low band dishes, 157 array beams from 
the coherent combination of dishes in the core are re- 
quired to equal the event rate of the incoherently com- 
bined total array. Coherent combination with dishes 
would only be optimal for hundreds of array beams, 
and even then it is likely that the processing power 
would be more effectively spent on AA-low. However, 
an incoherent commensal survey with low band dishes 
would be a cost-effective method to cover parameter 
space. 

5.2 Frequency and bandwidth effects 

The relative event rate per beam between signal com- 
bination modes described in the previous section is 
independent of frequency for searches of extragalac- 
tic populations and the first-order analysis of Galac- 
tic populations. The frequency-dependent effects in 
Equation 4 are common to all modes; the exception 
being processed FoV, however that cancels when the 
relative rate is calculated. 

For a given signal combination mode, the event 
rate per beam changes with frequency, as shown in 
Figure 3. For the extragalactic case (or more gener- 
ally, when the observed pulse width is much greater 
than the scatter broadening), the slope on the dish 
event rate is due to the spectral index of the source and 
changing FoV; for aperture arrays, the frequency de- 
pendence of Tays (due to sky noise) and Ac is also a fac- 
tor. Additionally, some pulsars display a turnover (a 
break in the spectrum where pulsar brightness is maxi- 
mum) around 100-200 MHz (Malofeev et al. 1994), and 
others above 1 GHz (Kijak et al. 2011). The event rate 
for such sources would decrease below the turnover fre- 
quency. For the Galactic case, the scatter broadening 
increasingly reduces the S/N for lower frequencies, in 
addition to the other factors described above. 

For SKAl, 7?.bcam-i is generally higher for AA-low 
than dishes, however the opposite is true for sources 
with low spectral indices and directions of larger scat- 
ter broadening. Also, lower frequencies have an in- 
creased memory cost for dedispersion. When low spec- 
tral indices, increased scatter broadening or low-frequency 
turnovers are factors and only a few AA station or ar- 
ray beams are available, the incoherent combination of 
the low band dishes gives a higher 7?.^. If SKAl dishes 
are equipped with PAFs, this would increase the total 
incoherent combination event rate of the dish array, as 
multiple station beams does for AA-low. 

The frequency and bandwidth effects on event rate 
per beam are interdependent. Section 4.4 shows that 
an optimal frequency range for searching for fast tran- 
sients with SKAl may be smaller than the full band 
of the receptor. As the frequency increases, the FoV 
reduces and the source luminosity is expected to de- 
crease for sources with pulsar-like emission character- 
istics. Processing more channels (hence bandwidth) 
increases TZ^, however this increase is less than Au^'''^. 

For some threshold beyond which extra bandwidth 
contributes little to the event rate, the processing for 
channels above this threshold could be more effectively 
used to form and search extra beams (increasing TZi, 
through FoV), trial more DMs or increase the detec- 
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tion S/N (sensitivity) through more optimal dedisper- 
sion techniques. The threshold depends on the costs 
of forming and searching the beams, which increases 
with bandwidth. Also, the number of channels (hence 
bandwidth) required to reach this threshold reduces for 
steeper spectrum sources (see Table 5). For searches of 
populations where the event rate is increasing with fre- 
quency, such as directions of large scatter broadening, 
the first channel is at the highest frequency and further 
contributing channels come from lower frequencies. 

Even if the search does not need the full band, 
recording it in a buffer is desirable to enable a more 
powerful analysis of detected transients with dedicated 
processing — this is especially important for relatively 
rare events. For flat spectrum sources, the S/N im- 
proves by a factor of V Au. The extra spectral infor- 
mation is also a useful analysis tool. For example, the 
multipath propagation which causes scatter broaden- 
ing would be more evident with a larger bandwidth, 
given the z/"* * relationship in Cordes & Lazio (2002), 
allowing us to distinguish between the intrinsic pulse 
width and a pulse that has been scatter broadened. 

5.3 Application to Galactic popula- 
tions 

A search for a Galactic population of fast transients 
is dependent on observing direction. The first-order 
analysis of scatter broadening shows that the event 
rate varies as a function of observation direction and 
intrinsic pulse width. The dependence of event rate 
on frequency, described previously, still applies. For 
an intrinsic pulse width of 1 ms, low band dishes give 
a higher event rate than AA-low for directions on or 
near the Galactic plane where the scatter broadening 
is large. Although the effectiveness of AA-low is re- 
duced at these directions, the increased luminosity of 
steep spectrum sources at lower frequencies somewhat 
counters this effect. Observations using the high band 
feed (1-2 GHz) on dishes are not modelled here, but 
in areas closer to the Galactic Centre, where scatter 
broadening is increased, it would achieve an event rate 
higher than the low band feed for spectrally shallow 
sources and sources with smaller intrinsic pulse widths. 

The event rate per beam curve for populations near 
the Galactic plane (Figure 3) is quite fiat in compar- 
ison to extragalactic populations, especially for shal- 
low spectrum sources. Because of this, additional pro- 
cessed bandwidth is more useful than for extragalac- 
tic populations, although it still cannot increase the 
event rate by more than the ideal case of Au^^*. For 
the threshold discussed in Section 4.4, more channels 
would contribute to the event rate. 

Besides scatter broadening, another limit on the 
event rate occurs if the telescope is sensitive enough 
to observe the population to the edge of the Galaxy 
(i.e. the sources are luminous enough to be observable 
to the edge of the Galaxy). In that case, the event rate 
for a sensitivity limited volume (as used in this paper) 
is not valid, because the limit is instead imposed by 
the boundary of the Galaxy. Macquart (2011) cap- 
tures both of these effects by introducing a direction 



dependent factor 5, such that 

7?.^ oc riprocS'^in eveuts s , (22) 

where < 5 < 3/2. A key result in Macquart (2011) 
is that FoV is more strongly preferred over sensitivity 
when scatter broadening or volume boundary limits 
increase S. Further modelling of the directional de- 
pendence of S is shown in that paper. 

The total number of events detected in a transients 
survey is proportional to TZ^t, where t is the total 
observation time. When S is large, the contribution 
from Smin is lessened. In this case, the total number 
of events detected can be more effectively increased 
by forming multiple beams (if available), incoherently 
combining the array or subarray signals (gaining a 
much larger FoV that coherent combination) or simply 
spending more time observing the sky. A low cost com- 
mensal survey using the incoherently combined array 
achieves these, making it even more effective for pop- 
ulations near the Galactic plane. 

To determine the number of events detectable for 
a Galactic population, simulations must account for 
all frequency-dependent effects on 7?.t,oam-i a func- 
tion of sky direction. Simulations for LOFAR pulsar 
searches (van Leeuwen & Stappers 2010) account for 
2sky (which modifies telescope sensitivity) as a func- 
tion of direction and scatter broadening as a function 
of direction and distance. For a population with a 
distribution of luminosities, the event rate is distance 
dependent. By numerically integrating along each line 
of sight, Macquart (2011) accounts for sources at each 
distance step which can no longer be detected below a 
certain sensitivity due to scatter broadening. Such an 
analysis is beyond the scope of this work, and requires 
assumptions to be made about the luminosity and spa- 
tial distribution of the objects, and their intrinsic pulse 
widths. 

5.4 Further work 

Maximising T^-bcam-i involves trade-offs between the 
receptor, signal combination mode, observing frequency 
and bandwidth for a given sky direction. This first- 
order analysis assumes 7?.cost-i oc TZ^^^^^-i. The next 
logical step is to describe the cost of forming and search- 
ing beams, to maximise 7?.(.ost-i . This problem is com- 
plex because beamforming and dedispersion process- 
ing costs are architecture specific. However, modelling 
these costs will give a better understanding of the opti- 
mal bandwidth and frequency to use for a given survey 
strategy and sky direction. Areas of work to further 
maximise 7?.(,ost-i are: 

• Trade-offs between station size and number of 
stations, for fixed AA-low collecting area. 

• The S/N gain versus processing cost from using 
a more accurate incoherent dedispersion method. 

• The S/N gain versus processing cost of using 
coherent dedispersion. Incoherent dedispersion 
has a lower S/N ratio than coherent dedisper- 
sion, but requires significantly less processing. 
To maintain the same detection rate as coherent 
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dedispersion, how many more beams need to be 
formed and searched? 

6 Conclusions 

Radio telescopes are inherently flexible, giving rise to 
a large design trade-off space, but any trade-offs must 
be made with the consideration of cost. This paper 
presents a new figure of merit, event rate per beam 
(^boam-i)i to measure the effectiveness of a survey 
strategy in detecting transient events in a volume of 
sky, while considering first-order costs — the ultimate 
goal being to calculate the event rate per unit cost 
(^cost-i)- The results show the complexities in deter- 
mining an optimal receptor, signal combination mode 
and frequency range for a given array. They highlight 
some important points for searching for fast transients 
with SKAl and more broadly with radio telescope ar- 
rays. 

The per beam event rate enables a frequency in- 
dependent analysis of the optimal signal combination 
mode and array configuration for a given receptor. 
This analysis is applied to SKAl, but regardless of the 
telescope used, incoherent combination always achieves 
a higher TZt^eam-^j making it preferable to coherent 
combination, subarraying and fly's eye in most cases. 
This is due to the event rate for surveying a volume 
of sky increasing faster with FoV than with sensitivity. 
The advantage of incoherent combination increases with 
a lower filling factor or the combination of an array 
with more elements. The exception to this result can 
occur once the number of dish or station beams formed 
reaches its maximum, due to physical and processing 
constraints. In this case, forming more array beams 
can result in coherent combination achieving a higher 
total event rate than incoherent combination. Fly's 
eye is unattractive because of the need to search many 
more beams, and its lack of localisation capability. 
However, if the search cost is low, incoherently com- 
bined three-element subarrays could usefully employ 
antennas unused by the primary (e.g. imaging) user 
observation. An example with SKAl would be split- 
ting 24 AA-low stations outside the core into 8 three- 
element subarrays, to achieve a total event rate similar 
to the incoherent combination of the total array. 

From Section 4.2 the per beam event rates for in- 
coherent and coherent combination can be compared 
using 

„3/2^3/4 

'^coh beam ^ ' ^'inc beam ^ 5 

where A'^o is the number of stations or dishes used for 
incoherent combination, rj is the fraction of A^o used 
for coherent combination and A^pix is the number of 
coherently formed array beams required to fill the FoV 
of the station or dish beam. This reduces to 

'*'cohbeam~^ ' *'inc beam^ 

for the ideal densely packed case, although the ratio is 
actually higher due to physical limitations. If the No 
antennas are divided into TVsa subarrays, 

'^sabeam~l -^'sa '^iiic bcam~ ■ 



Dependence on observing frequency and bandwidth 
needs to be considered when making event rate cal- 
culations, as does the observation direction and ex- 
pected spectrum of the source. The results in this 
paper show that for commensal observing, where the 
direction of observation is determined by the primary 
telescope user, a fast transients pipeline must dynam- 
ically adjust the search strategy to achieve the highest 
event rate. For a targeted survey, the choice of recep- 
tor needs careful consideration. 

For extragalactic searches with SKAl, the full avail- 
able bandwidth does not need to be searched; it is not 
an optimal use of the processing system. The contri- 
bution to event rate from processing additional band- 
width decreases as the signal frequency increases. Be- 
yond a threshold, the processing could be more effec- 
tively used to form and search more beams, trial more 
DMs or increase the detection S/N ratio. However, the 
full band is desirable for analysis of detected transients 
and for searches in directions where scatter broadening 
is such that the event rate is approximately constant 
with frequency. 

7 Recommendations for SKAl 

• SKA processing needs to provide flexible 
search modes 

The preferred receptor and signal combination 
mode depends on direction, especially for Galac- 
tic populations. Incoherent and coherent combi- 
nation modes are both effective with SKAl and 
depend on the array filling factor, signal process- 
ing costs and the spectrum of the source popu- 
lation. 

For SKAl, coherent combination of the AA-low 
core requires approximately five array beams to 
equal the event rate of a single beam of the in- 
coherently combined total array. For low band 
dishes, this number is 157, making the coherent 
combination of dishes unattractive unless more 
than this many array beams can be formed. 
SKAl AA-low is effective for extragalactic tran- 
sients searches (assuming dedispersion costs are 
not too high), especially for steep spectrum sources 
(see Figure 3). For fiat spectrum sources, the 
advantage of aperture arrays disappears. Low 
band dishes give a higher event rate for direc- 
tions near the Galactic plane, where scatter broad- 
ening is larger. 

• The re-use of signal processing required 
for SKA imaging modes enables low cost 
fast transients searches 

Commensal surveys using the incoherent com- 
bination of dish signals or AA station beams 
are low cost options for searching for fast tran- 
sients with SKAl. The incremental cost of im- 
plementing such a search is small because it uses 
beams formed for the primary user observation. 
A commensal survey using both receptors is a 
simple method to increase the total number of 
events detected, by increasing the total observa- 
tion time. The survey effectiveness, compared to 
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the coherent combination of the array, increases 
for populations closer to the Galactic plane. 

Access to the dish and station beam data should 
be via a spigot with a well defined interface, 
to enable the implementation of flexible search 
modes. The processing for commensal surveys 
could be implemented internally or with user- 
provided processing as it becomes affordable or 
available. 

• Requirements for fcist transient searches 
with SKAl 

Until the processing costs are further explored, 
the base requirements for SKAl are: 

— Availability of incoherent and coherent com- 
bination modes for AA-low and low band 

dishes. 

— Processing for low cost commensal survey 
modes; or provision for access to the dish 

and AA station beam data via spigots. 

— Voltage (coherent) buffering capability of 
the full band; of order tens of seconds for 
dish frequencies and possibly minutes for 
lower frequencies. 

For extragalactic searches, processing the full 

available bandwidth is not required. Bandwidths 
of 50- 100 MHz arc sufScicnt on the basis of the 
simplified investigation undertaken in this pa- 
per; a more detailed study of the trade-offs could 
be made. However, buffering the full band is de- 
sirable for analysis of detected transients. 
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The following tables show how sensitivity, FoV, beam- 
former processing cost and the number of data streams 
to be searched scale for different signal combination 
modes. The scaling equations assume that the polari- 
sations are summed prior to searching and all elements 
are of equal diameter and sensitivity. 

A. 2.1 Incoherent combination 

Incoherent combination of the element signals requires 
the signal from each element to be detected, a geomet- 
ric delay applied and the signals summed. Figure 5 
shows the steps to incoherently combine signals and 
Table 6 shows some performance attributes. 



RF/IF -J 



Averager 



I- I 



Figure 5: Incoherent combination. For aperture 
arrays, station beamforming takes place prior to 
the y2 block. 



Appendices 

A Signal combination techniques 



Table 6: Incoherent combination attributes 



The performance attributes of a radio telescope array 
depend on how the signals from the array elements are 
combined. This in turn affects the detection rate for 
fast transients. The array elements may be an antenna 
(such as a dish or dipole) or a phased group of antennas 
(stations). We term the single antenna primary beam 
or the phased station beam as the element beam, with 
FoV Qq. The signals from these elements may then be 
combined incoherently or coherently as discussed later 
in this section. 

A.l A A station beamforming 

An aperture array station will have thousands of indi- 
vidual antennas. To reduce the data rate from the 
station and the downstream signal processing load, 
the antennas will be phased into a station beam, as 
described in Zarb Adami et al. (2010). Multiple sta- 
tions beams can be formed by applying the appropriate 
phase shift. The beamformer processing cost scales as 
A'^a/st A'b-o Ai^ operations per second, where Ns,/st is 
the number of antenna elements per station, A^b-o is 
the number of station beams formed and Af is the 
bandwidth (Cordes 2009). 

A. 2 Modes of beamforming for search- 
ing 

For antenna, phased array feed or station beams point- 
ing at the same location on the sky, the signals de- 
tected may be combined incoherently or coherently. 



Attribute 


Scaling 


Comment 


Sensitivity 




Nq elements, each 






with effective area 






Aco 


FoV {rtproc) 




Nh-o beams 






formed per station. 






each with FoV ilo'' 


Processing 






Data streams 







^ rig = I (^S^) ' where /Co is the dish illumina- 
tion or station beam taper and Dq is the diameter 
of the dish or station in metres. 



As long as the appropriate geometric delay is ap- 
plied to the signal at each element, incoherently com- 
bined elements do not need to be located close to- 
gether. For the SKA, this means that while the core is 
being used for low angular resolution experiments, the 
mid and long baselines could be used for fast transient 
searches. Note that incoherent combination cannot ac- 
count for the geometric delays within the beam, but 
away from the beam centre. 

We assume that digitising, channelising and sta- 
tion beamforming are existing telescope functions and 
that these do not factor into the additional process- 
ing cost of incoherently combining the station beams. 
Square-law detection involves squaring and summing 
the real and imaginary components of each channel of 
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each beam from each station. The averager integrates 
the power samples for each beam, and then equivalent 
beams from different elements are summed together, 
resulting in a total of Nh-o incoherently combined 
data streams. The integration of the power samples 
greatly reduces the data rate of the incoherent beams, 
although it comes at a cost of a lower time resolution. 




Averager — - 



Figure 6: Coherent combination. For aperture ar- 
rays, station beamforming takes place prior to the 
X; block. 



A. 2. 2 Independently pointed subarrays, in- 
coherently combined 

A further increase in FoV can be achieved by pointing 
subarrays of elements in different directions and inco- 
herently combining the signals from elements in the 
subarray. Table 7 shows performance attributes for 
this mode. The processing cost for incoherently com- 
bining subarrays is approximately the same as it is for 
incoherently combining all stations because the same 
number of beams need to be square-law detected and 
summed across stations; the difference is that separate 
sums need to be maintained for each subarray. The 
approximate beamforming operations cost is therefore 
independent of the number of incoherently combined 
subarrays (A^'sa). 



Table 8: Coherent combination attributes. 



Attribute 


Scaling 


Comment 


Sensitivity 




Coherent sum of 






Nq elements 


FoV (riproc) 


AT O ^ 


A'b-arr array beams 






formed 


Processing 




Data streams 






Data streams 


iVb-oA^pix 


A'pix pixels to fill 


to pixelise full 




rig with array 


FoV 




beams (see 



Equation 5) 



Table 7: Incoherently combined subarray at- 
tributes. 



2 la 



TT / C/Carr A 
4 \^!^-Darr^ 



where /Carr is the array beam 



taper and I?aiT is the longest baseline in the array. 



Attribute 


Scaling 


Comment 


Sensitivity 




No/sa elements per 


(A) 


subarray 


FoV (riproc) 




Nsa independently 






pointed subarrays 


Processing 






Data streams 


iVsaiVb-O 





The smallest subarray size is one element, meaning 
each element is pointing to a unique patch of sky — 
this is termed fly's eye. In reality, the minimum num- 
ber of elements in a subarray would be 3, to allow for 
a triggered buffer to be used to localise any detected 
signal. 



A. 2. 4 Correlation beamforming — 'fast imag- 
ing' 

An alternative method to pixelise the full FoV is to 
create images from the correlator at a high time resolu- 
tion. The highest time resolution available is specified 
by the correlation integration period, which must be 
short enough to prevent smearing. Cordes (2009) com- 
pares the computational cost of this method with the 
coherent sum, and finds that the cost-effectiveness de- 
pends on configuration and system temperature. How- 
ever the major advantage is that the correlator hard- 
ware can be re-used. It should be noted that this 
means that the data must be able to be dumped from 
the correlator at these high rates which may place ad- 
ditional requirements on the correlator architecture. 
There is also a computational cost for gridding and 
imaging at high time resolution. 



A. 2. 3 Coherent combination — array beam- 
forming 

For coherent combination, each array beam is formed 
by the weighted sum of A'^o element beams pointing in 
the same direction. Geometric delays are applied to 
the signals from the elements, which are then summed 
and detected. Like AA stations, multiple beams can 
be formed. Figure 6 shows the steps for coherent array 
beamforming, and Table 8 shows the performance at- 
tributes. Cordes (2009) discusses array beamforming 
in the context of the SKA. 



Table 9: Correlation beamforming attributes. 



Attribute 


Scaling 


Comment 


Sensitivity 




For large A'^o- 


(A) 






FoV (riproc) 






Processing 


1 


Use existing 






correlator hardware 


Data streams 







18 



Publications of the Astronomical Society of Australia 



B Event rate for a broadened 
pulse in a volume of sky 

The detected event rate is calculated for observable 
space. We assume extragalactic sources of intrinsic 
luminosity d Jy pc'^ are homogeneously distributed in 
a sphere of volume V with a nominal event rate density 
of peventss"^ pc~^. The event rate for this sphere is 
given by pV: 



TP _ ^ 

''^sphere ^3 



events s 



(B.l) 



where D is the radius of the sphere in pc. From Deneva 
et al. (2009), a source can be detected out to a maxi- 
mum distance 



^max — 



Li 



(B.2) 



where Si is the intrinsic flux density of the source. 

The intrinsic flux density of a pulse differs from 
the observed flux density due to pulse broadening (or 
smearing) effects of the interstellar medium and of the 
detection system itself. However, if no energy is lost 
due to attenuation, then pulse "area" is conserved such 
that SiWi — SW, where Wi is the intrinsic width of 
the pulse, S is the observed flux density and W is the 
width of the broadened pulse. If the post-detection in- 
tegration time r is equal to W, a telescope with mini- 
mum detectable flux density Sinin,T^w can detect the 
broadened pulse to a maximum distance 



^max — 



WiCi 



4ttWS, 



mill. T — W 



(B.3) 



A more general relationship which uses an integration 
time of T = Wi is 



Ci 



(B.4) 



Pulse broadening factors are discussed further in 
Appendix C. 

The nominal extragalactic population is observ- 
able out to Dmax for the fraction of the sky observed 
(fiproc/f2sky) aud the detected event rate is 



47r ^^proc 



3 ^ 



^sky 



— 3 



( 



V 47r5'min 



(B.5) 



f^proc, Smin and W are functions of how the signals 
are combined and processed by the telescope system. 
•Smin = crAS'rms, where a is the required S/N ratio 
and A5rms is the rms variation in flux density for a 
randomly polarised source, obtained by applying the 
telescope gain to the radiometer equation: 



5'min ^AiS'rms 



(B.6) 



where Tsys is the system temperature, A^ is the to- 
tal effective area used in the survey, Ai/ is the pro- 
cessed bandwidth and r is the post-detection integra- 
tion time, which also defines the time resolution. 



C Pulse broadening and cor- 
rection (dedispersion) 

Cordes & McLaughlin (2003) model the broadening 
of a delta function pulse due to propagation through 
the interstellar medium and signal processing response 
times using the following approximation: 



At = 



(C.l) 



where AtoM represents the broadening due to disper- 
sion smearing; AtaoM is due to the error in the DM, 
SDM, used by the system's dedispersion signal pro- 
cessing; AtAv is the system's filter response time; and 
Td is due to the multipath propagation effects of the 
medium. Expanding on this, the broadening of a pulse 
of intrinsic width Wi can be modelled as 



W = 



/W^ + Atl^ 



The filter response time, AIav, is approximately 
equal to Ai'~^ , where Au is the bandwidth of the fil- 
tered signal. It is important to note that for a fully 
coherent transient detection system where the element 
beams are coherently combined, coherently dedispersed 
and searched, Au represents the full signal bandwidth; 
whereas for an incoherent transient detection system in 
which the signal is channelised, detected and searched, 
Av represents the much smaller channel bandwidth. 
Consequently, the filter response component of pulse 
broadening, AtAv, is significantly higher for incoherent 
transient detection systems. 

Furthermore, while coherent dedispersion techniques 
can completely correct for dispersion smearing (given 
that the DM is known), incoherent dedispersion tech- 
niques can only correct for dispersion between the filter- 
bank channels; they cannot correct for dispersion within 
the channels. Intra-channel dispersion smearing can 
be reduced by choosing smaller channel bandwidths, 
but at the expense of larger filter response times. The 
optimum channel bandwidth for incoherent dedisper- 
sion occurs where the dispersion smearing within each 
channel equals the filter response time (Hankins & 
Rickett 1975; Cordes & McLaughhn 2003). This leads 
to a minimum pulse width after incoherent dedisper- 
sion of 



Wine = s/Wi + 2(AtDMmin)2 "f Atj^^ + tJ , (C.3) 

where AtDMmin = VS. 3 x 10^^ DMiy-'"^ . It should be 
noted that this optimum cannot be realised for all 
DMs, because it expects the channel bandwidth to be 
a function of the DM. 

For coherent dedispersion the AtuM term is com- 
pletely removed: 



Wcoh = ^/W^ + Atl^ + AtioM + ri, (C.4) 

where (as noted above) Av is the full processed band- 
width and AtAi/ ~ Av~^. 
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D Frequency dependence of low 
frequency aperture arrays 

For illustrative purposes, we plot the event rate TZv and 
a breakdown of its frequency-dependent components. 
Figure 7 shows TZi, for source luminosities with spectral 
indices ^ = —1.6 and 0, over a frequency range 70- 
450 MHz, at 1 MHz steps with processed bandwidth 
/\v = \ MHz and normalised to = 1 at 70 MHz. 
The slope is steep: at 160 MHz, the event rate is 10% 
of the event rate at 70 MHz. At 450 MHz, the event 
rate is 0.016% of the 70 MHz event rate. A further 
breakdown of Smin^ is shown in Figure 8. 



E Event rate as a function of 
frequency 

The event rate is determined numerically by calculat- 
ing the S/N ratio (SNR) for each of A^ch frequency 
channels of width Ai^ch, weighting it by w, the amount 
of sky seen with that SNR and then taking the root of 
the sum of the squares: 



\ 



^{Cw^SNR,f, (E.l) 



where C is a constant. Given SNR = 

4-kS ' 

2 

Setting w = Hp^oc and C = {^pf''^ so TZl^^ = Cxwx 
SNR, we can say that 



J- / *— 1 I / proc,z z 



E 



(E.3) 



for a processed bandwidth of Ai/ = AchAi^ch. Where 
there is no frequency dependence, the event rate for 
total bandwidth and unit bandwidth Ai^ch = 1 
becomes 

3 

7^A. = JpA^^'/'^fiproc ( ^ ) ' (E.4) 



as expected. 
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Figure 8: S'min„ and breakdown of the frequency-dependent components comprising S'min„, normalised to 
the maximum value of each. 



